THE DIFFUSIVE COMPETITION MODEL WITH A FREE BOUNDARY: 
INVASION OF A SUPERIOR OR INFERIOR COMPETITOR* 



YIHONG DUt AND ZHIGUI LIN* 

Abstract. In this paper we consider the diffusive competition model consisting of an invasive 
species with density u and a native species with density v, in a radially symmetric setting with 
free boundary. We assume that v undergoes diffusion and growth in R^, and u exists initially in 
a ball {r < /i(0)}, but invades into the environment with spreading front {r = h{t)}, with h(t) 
evolving according to the free boundary condition h' [i) = —fmr{t,h{t)), where ^ > is a given 
constant and u{t,h{t)) = 0. Thus the population range of u is the expanding ball {r < h{t)}, 
while that for v is R^. In the case that w is a superior competitor (determined by the reaction 
terms), we show that a spreading- vanishing dichotomy holds, namely, as t — oo, either h{t) — s- oo 
and {u,v) — >■ (u*,0), or Vmit^ooh{t) < oo and {u,v) — >■ (0,i;*), where (u*,0) and {0,v*) are the 
semitrivial steady-states of the system. Moreover, when spreading of u happens, some rough 
estimates of the spreading speed are also given. When u is an inferior competitor, we show that 
(u, v) — >■ (0, V*) as t ^ oo. 



1. Introduction 

In this paper we study the behavior of the solution (u{t,r),v{t,r),h(t)) to the following 
reaction-diffusion problem with radial symmetry, 

ut — diAu = u{ai — biu — civ), t > 0, < r < h{t), 

vt — d2Av = v{a2 — — C2v), t > 0, < r < oo, 

Ur{t, 0) = Vrit, 0) = 0, u{t, r) = 0, t > 0, h{t) < r < oo, 

h'{t) = -flUr{t,h{t)), t>0, 

hiO) = ho, n(0, r) = no(r), < r < /iq, 

w(0,r) = VQ^r), < r < oo, 



(1.1) 



where An = Urr + ^^^y^Ur, r = h{t) is the moving boundary to be determined, /iq, /i, di, ai, bi 
and Ci{i = 1,2) are given positive constants, and the initial functions no and vq satisfy 

. . / Mo e C2([0, ho]), u'q{0) = no(/io) = and no > in [0, ho), 

^ \ vo€ C2([0, +oo)) n L°°(0, +oo), v'oiO) = and no > in [0, +oo). 

Ecologically, this problem describes the dynamical process of a new competitor invading into 
the habitat of a native species. The first species (n), which exists initially in the ball {r < ho}, 
disperses through random diffusion over an expanding ball {r < h{t)}, whose boundary {r = h(t)} 
is the invading front, and evolves according to the free boundary condition h'{t) = —^Ur{t, h{t)), 
where // is a given positive constant. The second species (n) is native, which undergoes diffusion 
and growth in the entire available habitat (assumed to be MJ^ here). The constants di and d2 
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(1.3) 



are the diffusion rates of u and v, respectively, oi and 02 are the intrinsic growth rates, hi and 
C2 are the intraspecific and ci and 62 the interspecific competition rates. 

In the absence of a native species, namely v = 0, the system reduces to the following diffusive 
logistic problem, 

Ut - diA = u{ai - biu), t > 0, < r < h{t), 
Ur{t,0) = 0, u{t,h{t)) =0, t>0, 
h'{t) = -flUr{t,h{t)), t>0, 
hiO) = ho, uiO, r) = uoir), < r < /iq, 

which has been treated in [4], extending the one dimensional case first studied in [5]. The 
behavior of (|1.3p is characterized by a spreading- vanishing dichotomy, namely, as t — )• 00, one of 
the following alternatives occurs: 

• Spreading: h{t) — )• 00 and u{t,r) — )• ai/61, or 

• Vanishing: h{t) — )• h^o < 00 and u{t, r) — )• 0. 

Moreover, when spreading occurs, it is shown that h{t)/t — )• fco G (0, 2^/ aidi) as t — )• 00, and /cq 
is called the asymptotic spreading speed of u. Further discussions of and a deduction of the 
free boundary condition based on ecological assumptions can be found in [1]. 

In this paper, we will examine the case that u invades into an environment where a native 
competitor already exists. This is a much more complicated situation, and we will only consider 
(jl.ip under certain restrictions on the parameters, to be specified below. 

Problem p.ip is a variation of the diffusive Lotka-Volterra competition model, which is often 
considered over a bounded spatial domain with suitable boundary conditions or considered over 
the entire space ([HI [l3]). For example, the dynamical behavior of the following bounded 
domain problem 

Ut — diAu = n(ai — hiu — civ), {t, x) E (0, 00) x il, 

Vt — d2Av = v{a2 — b2U — C2v), {t, x) G (0, 00) x il. 



(1.4) 



n(0, x) = no(x) > 0, f(0, x) = vo{x) > 0, x £ 



is well known, where is a bounded smooth domain of with N > 1, is the outward 
unit normal vector on dfl. This model describes the situation that two competitors evolve in a 
closed habitat with no flux across the boundary d^l. Therefore their competitive strengths 
are completely determined by the coefficients (oj, hi, Ci, di) in the system, i = 1,2. 

Problem ()1.4p admits the trivial steady state Rq = (0, 0) and semi-trivial steady-states i?i = 
(ai/61, 0) and i?2 = (0, 02/02). Moreover, if 61/62 > «i/«2 > ci/c2 or 61/62 < ai/«2 < ci/c2, 
the problem has a unique constant positive steady-state 

j^* _ f Q1C2 - a2Ci 0261 - 0162 ' 

V 61C2 - 62C1 ' 61C2 - 62C1 

These are all the nonnegative constant steady-states. There may also exist non-constant positive 
steady-states, but they are all linearly unstable when is convex ([9]). For the constant equilibria, 
their roles are summarized below (see, for example, |13j page 666): 

(1) Rq is always unstable; 

(2) when 61/62 > 0^1/02 > Ci/c2, R* is globally asymptotically stable; 

(3) when a\ja2 > max{6i/62, C1/C2}, Ri is globally asymptotically stable; 

(4) when 01/02 < min{ci/c2, 61/62}, R2 is globally asymptotically stable; 

(5) when 61/62 < 01/02 < Ci/c2) -^i and R2 are locally asymptotically stable, and R* is unstable. 
In case (2), the competitors co-exist in the long run, and it is often referred to as the weak com- 
petition case, where no competitor wins on loses in the competition. In case (3), the competitor 
u wipes V out in the long run and wins the competition; so we will call u the superior competitor 
and V the inferior competitor. Analogously u is the inferior competitor and v is the superior 
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competitor in case (4). Case (5) is the strong competition case, and the long-time dynamics of 
(II. 4p is usually complicated and difficult to determine. 

We will only consider cases (3) and (4) for (jl.ip . We will show that in case (3), similar to 
(jl.3p . a spreading- vanishing dichotomy holds for (|l.ip . namely as t — )• oo, either h{t) — )• oo and 
{u,v) — )• Ri (spreading for u), or h{t) — )• hoo < oo and {u,v) — )• R2 (vanishing for u). Clearly this 
is strikingly different to the long-time behavior of (jl.4p . However, in case (4), we show that as 
t — )• 00, {u,v) — )• R2, so the dynamical behavior is similar to that of (jl.4p in case (4). 

For the entire space problem 

, . j ut — diAu = u{ai — biu — civ), {t, x) G (0, 00) x M^, 

\ vt-d2Av = v{a2-b2U-C2v), (t, x) G (0, 00) x M^, 

extensive work has been done concerning the existence of traveling wave solutions in space di- 
mension N = 1. For example, for case (3), it is shown in [8] that there exists c* > such that 
for each c > c*, (11. 5p with = 1 has a solution of the form 

(1.6) (n(x, t),v{x, t)) = {U{x - ct),V{x - ct)) 
satisfying 

(1.7) U' <0,V' > 0, {Ui-oo),V{-oo)) = Ri, {Ui+oo),V{+oo)) = R2] 

there is no such solution when c < c*. The general long-time behavior of the Cauchy problem of 
(II. Sp . however, is still poorly understood (see Remark 13.21 below for a partial result). 

We end the introduction by mentioning some related research. In [Tl], a predator-prey model 
in one space dimension was considered, where the available habitat is assumed to be a bounded 
interval [0,/], and no-flux boundary conditions are assumed for both species, except for the 
predator sX x = I. It is assumed that the predator satisfies a free boundary condition as in (II. ip . 
before the free boundary x = h{t) reaches x = I, and a no-flux boundary condition at x = I 
is satisfied by the predator after the free boundary has reached x = I. It is shown in [11] that 
the free boundary always reaches I in finite time, and hence the long-time dynamical behavior 
of the free boundary problem is the same as the fixed boundary problem. After the first version 
of this paper was completed, we have learned several more closely related research. In [7], the 
week competition case was considered in one space dimension, but in their model, both species 
share the same free boundary. Such a free boundary setting was also used in jl5j for the Lotka- 
Volterra predator-prey system in one space dimension. In |16] . the Lotka-Volterra predator-prey 
model was considered in one space dimension, where similar to (jl.ip . one species (the predator) 
is subject to free boundary conditions, and the other is considered over the entire M^. 

The rest of this paper is organized in the following way. In section 2, we prove some general 
existence and uniqueness results, which implies in particular that (jl.ip has a unique solution 
defined for all t > 0. Moreover, some rough a priori estimates are given, as well as a rather 
general comparison result. These results are useful here and possibly elsewhere. In section 3, 
we investigate the case that u is an inferior competitor, namely the coefficients fall into case (4). 
Sections 4 and 5 are devoted to the case that n is a superior competitor. A spreading-vanishing 
dichotomy is established in section 4, and a sharp criterion to distinguish the dichotomy is also 
given there. In section 5, some rough estimates for the spreading speed is given for the case that 
spreading of u happens. 

2. Preliminary Results 

In this section, we first prove a local existence and uniqueness result for a general free boundary 
problem, and then we obtain global existence results, which imply that the solution to (II. ip exists 
for all time t G (0, 00). Lastly, we obtain some comparison results, which will be used in the other 
sections. 
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Consider the following general free boundary problem: 

ut - diAu = f{u,v), t > 0, <r < h{t), 

vt — d2Av = g{u,v), t > 0, < r < oo, 

Ur{t, 0) = Vr{t, 0) = 0, u{t, r) = 0, t > 0, h{t) < r < oo, 

h'{t) = -^lUr{t,h{t)), t>0, 

h{0) = ho,uiO,r)=uo{r), < r < /iq, 

f(0,r) = VQ{r), < r < oo, 

where /(O, v) = g{u, 0) = for any u,v gR, and uq, vq are as in (jl.ip . 

Theorem 2.1. Assume that f and g are locally Lipschitz continuous in M^. For any given 
{uq,vq) satisfying ()1.2p and any a £ (0,1), there is a T > such that problem (j2.ip admits a 
unique bounded solution 

moreover, 

(2.2) ||n||p(l + Q)/2,l + Q(£)^-) + ||t'||c'(l + Q)/2,l + c«(-2)oo) + ||/i|lcl + a/2([0,r]) — ^' 

w/iere L>t = {{t,r) e : t E [0,r],r G [0,/i(t)]}, = {(t,r) e : t e [0,T],r G [0,+oo)}, 
C and T only depend on Hq, a, ||^io||c2([o,/io])' Il^o||c72([o,oo)) '^'^^ ^he local Lipschitz coefficients of 
f,9- 

Proof: The proof is similar to that in [5j and [4J for the scalar problem, with some modifica- 
tions. We sketch the details here for completeness. First we straighten the free boundary as in 
[2]. Let ({s) be a function in (7'^[0, oo) satisfying 

C(s) = 1 if |s - ho\ < ^, C(s) = if |s - ho\ > ^, |C'(s)| < ^ for all s. 



ha 



Consider the transformation 

{t,y) ^ (t,x), where x = y + C{\y\){h{t) - /iq)^, y G i?^, 

\y 



which leads to the transformation 

(t, s) —?- {t, r), with r = s + C{s){h{t) - ho), < s < oo. 

As long as 

\Ht)-ho\<^, 

o 

the above transformation x — )• y is a diffeomorphism from R''^ onto and the transformation 
s — 7- r is also a diffeomorphism from [0, +oo) onto [0, +oo). Moreover, it changes the unknown 
free boundary |x| = h{t) to the fixed sphere \y\ = ho. Now, direct calculations show that 



Let us also denote 



dr 1 + C{s){h{t) - ho) 
d^s _ C'{s){h{t) - hp) 
dr^ [l + C'is){h{t)-ho)]^ 

1 ds _ C(g) 
'W{t)di ~ l + C'(s)(/i(t)-/io) 

{N -l)^/A 



Bihit),s), 
C{h{t),s). 



s + C{s)ih{t) - ho) 



D{h(t),s). 
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If we set 

u{t, r) = u{t, s + C{s){h{t) - ho)) = w{t, s), 
v{t, r) = v{t, s + C{s)iMt) - ho)) = z{t, s), 

then the free boundary problem ()2.ip becomes 



(2.3) 



wt — AdiWss — {Bdi + h'C + Ddi)ws = f{w, z) 
zt - Ad2Zss - {Bd2 + h'C + Dd2)zs = g{w,z), 
Ws{t,0) = Zs{t,0) = w{t,r) = 0, 
h'{t) = -fiWs{t,ho), 
h{0) = ho, w{0, s) = wois) := uo(s), 
z(0,s) = zo{s) := vo{s), 



t>0, < s <ho, 
t > 0, < s < oo, 
t > 0, Hq < s < oo, 
t > 0, 

0<s<ho, 
< s < oo, 



where A = A{h{t),s), B = B{h{t), s), C = C{h{t),s) and D = D{h{t), s). 
We denote h* = —^u'o{hQ), wq{s) = for s > /iq and for < T < -^(^fh*)^ 

At = [0, T] X [0, ho] , A?? = [0, T] X [0, +oo) , 

HT = [he Ci[0,r] : /i(O) = ho, h'{0) = h\ \\h' - h*\\ci[o,T]) < l}, 
Wr = |u> G C{Af) : w{t, s) = for s > /iq, < t < T, 

w{0,s) = wo{s) for < s < Hq, \\w - -woUl^CA^?) = ll""^ - wo\\c{At) ^ l}i 
Zt = [z£ C{Af) : z{0,s) = zo{s), |k - ^o||l-{A-) ^ l}- 
It is not difficult to see that Tt '■= Wt x Zt x Ht is a complete metric space with the metric 

V{{wi,Zi,hi), {w2,Z2,h2)) = \\wi - W2\\c{At) + 11^1 - ^2||l°°(A??) + ll^'l " ^2 II C([0,T]) • 

Let us observe that for /ii, /i2 G Ht, due to /ii(0) = /i2(0) = ho, we have 

(2.4) \\hi - h2\\ci[0,T]) < T\\h[ - h'2\\ci[0,T])- 

Next, we shall prove the existence and uniqueness result by using the contraction mapping the- 
orem. Since / and g are locally Lipschitz continuous, there exists an L* depending on ||uo||c([o,/io]) 



and ||uo||l°°{[o,+oo)) such that 

\f{w,z)\ = \f{w,z) - f{0,z)\ < L*\w\ < L*(||no||c[o,M + 1)' ("''^) ^ ^t x Zt, 
\g{w, z)\ = \g{w, z) - g{w, 0)| < L*\z\ < L*(||'t;o|Uco[o,oo) + 1), {w, z) e Wt x Zt- 

By standard theory and the Sobolev imbedding theorem [10], for any {w,z,h) G Tt, the 
following initial boundary value problem 

Wt - AdiWss - {Bdi + h'C + Ddi)ws = f{w, z), t > 0, < s < ho, 

Zt — Ad2Zss — {Bd2 + h'C + Dd2)zs = g{w, z), t > 0, < s < oo, 

(2.5) I Ws{t, 0) = Zs{t, 0) = w{t, r) = 0, t>0, ho<r <oo, 

w{0, s) = wo{s) := uo{s), < s < ho, 

z{0, s) = zo{s) := vo{s), < s < oo 

admits a unique bounded solution {w,z) G C^^+^^/^.i+a^^^) ^ Q(i+a)/2,i+a^^^^ ^^^^ 
(2-6) ||w||c{i+Q)/2,i+a(Aj,) ^ Ci, 



(2.7) 



l^llc(i+")/2.i+'»(A3?) ^ ^li 
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where Ci is a constant depending on a,hQ,L*, Huollc^fc/io] \\'^o\\c'^[o,+oo)- The estimate ([27 
comes from the interior estimate. For any m > 0, by classical parabolic regularity theory [10] . 
one then have the estimate 

Pllc(i+")/2'i+"{[0,T]x(m,m+l)) ^ C'(a) || Z || vi/2,i,P([0,T] X [m,m+l]) ^ C{a,p,L* ,T) 

for some large p > 1 and a G (0, 1). 

Now, we define h{t){> 0) by the fourth equation in (j2.3p : 

(2.8) h{t) = ho-fi [ Ws{T,ho)dT, 

Jo 

which infers h'{t) = -fiWsit,ho), h{0) = ho and h'{0) = -/i'u^(/io) = h* . Hence h' £ C°/2([o,r]) 
with 

(2.9) ll^'llc(i+")/2{[o,T]) ^ := ^Ci. 
Now we define the map 

T: C{A^) X C{A^) x C^[0,r] 

by F{w{t, s),z{t, s); /i(t)) = {w{t, s), z{t, s); /i(t)). It's easy to see that {w{t, s), z{t, s); h{t)) £ Tt 
is a fixed point of J- if and only if it solves ()2.3p . 
The estimates in ([221), (1121) and ^M) yield 

11^' - h*\\c{[o,T]) < ll^'llc'»/2([0,T])^"^^ ^ ^CiT"/^, 

- W'ollcCA^) < IIW' - W;olb{i+-)/2,0(A^)T(l+")/2 < CiT(1+")/2, 

115 r II ^11; ^11 T^(l+a)/2 ^ n T^(l+")/2 

11^ — ^0||l°°{A5?) S IP - 2:ollc{i+<»)/2,0(A;^)-t S L-l^ ' 

Therefore if we take T < min{(^C7i)-2/"^ 2/{i+a)^^ ^^^^ ^ 

maps Tt into itself. 

Now we prove that for T > sufficiently small, is a contraction mapping on Tt- Indeed, let 
{wi,Zi,hi) G Ft (i = 1,2) and denote {wi,Zi,hi) = T{wi, Zi,hi). Then it follows from (j2.6p . (j271 
and (USD that 



ll^^i|lc(i+a)/2.i+a(AT) — C*!, Pi|lc{i+a)/2,i+a(A5?) — C'l, ||^illc"/2{[0,T]) — C'2- 

Setting W = wi — 11)2, we find that VF(t, s) satisfies 

Wt - A{h2, s)diWss - {B{h2, s)di + h'^C{h2, s) + D{h2, s)d{)Ws 

= [A{hi,s) - A{h2, s)]diwi^ss + [B{hi,s) - B{h2, s) + D{hi,s) - L»(/i2, s)]diwi,s 
+ [h[C{hi,s) - /i2C(/i2, s)]wi^s + f{wi,zi) - f{w2, Z2), t>0, 0<s <ho, 

dW 

— (t,0) = 0, Wit,ho) = 0, t>0, 
W{0,s) = 0, 0<s< /lo- 
using the estimates for parabolic equations and Sobolev's imbedding theorem, we obtain 

< C'adki - ?^2||c(At) + Iki - ^2||l«=(A5?) + 11^1 - ^2||ci([0,T])), 

where C3 depends on Ci, C2, the local Lipschitz coefficients of /, g and the functions A, B and C 
in the definition of the transformation (t, s) — )• {t,r). Similarly, we have 

I2ll\ Ikl - ^2|lc{i+'»)/2,i+<»{A5?) 

< C4(||wi - 'U;2||c{At) + 11^1 - ^2||l°o(A5?) + 11^1 - ^2||ci([0,T])), 
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where C4 depends on Ci,C2, the local Lipschitz coefficients of f,g and the functions A,B and 
C. Taking the difference of the equations for hi and /i2 results in 

(2.12) \\h'i - /i2llc"/2([0,T]) ^ /^(ll^l.s ~ W'2,s|lc<»/2,0(Ay))- 

Combining ([23]), (j2J0]) . (pTT]) and ([232]) . and assuming T < 1 we obtain 

ll'fi'l — '"^2|lc(i+a)/2.i+"(AT) + 11^1 ~ ^2|lc(i+a)/2,i+Q(A;^) + im ~ ^2llc"/2([0,T]) 
< Csdluii - U'2||c7(Aj,) + ll^i - 2:2||l°°(A5?) + ll^'l - ^2llc[0,T]), 

with C5 depending on 0^,04 and ^. Hence for 

{l + h*) 
we have 

ll'W'l - W2\\c{At) + Pl ~ ^2||L°°(A5f) + 11^1 - ^2llc([0,T]) 

< r'^-^+")/^(||7Z;i - '(i)2|lc(i+«)/2.i+"{AT) + 11^1 ~ ^2|lc(i+")/2>i+«(A5?)) + T'^'^^ll^i - ^2llc«/2([0,T]) 

< CsT^/^dlu;! - 'ii;2||c(AT) + Iki - ^2||L°°(A5f) + \\h'i - /i2l|(7([o,r])) 

< ^(||u;i - U'2||c(At) + Pi - ^2||l-'(A5?) + - /i2llc([0,T]))- 

This shows that for this T, is a contraction mapping in Tx- It follows from the contraction map- 
ping theorem that has a unique fixed point {w, z, h) in Tt- In other words, {w{t, s), z{t, s); h{t)) 
is the solution of the problem (|2.3p and therefore (u(t, r),v{t, r); /i(t)) is the solution of the prob- 
lem (|2.ip . Moreover, by using the Schauder estimates, we have additional regularity of the 
solution, h{t) G Ci+"/2[0,r], u G CiW2,2+a((-y) and v G C7i+"/2.2+"((o, T] x (0,+oo)). Thus 
{u{t, s),v{t, s); h(t)) is the classical solution of the problem ()2.ip . where Gt := {(i)^') G M2 : t G 
(0,T],r G (0,/i(t))}. □ 

Remark 2.2. By a bounded solution (u, f , /i) we mean that there exists Mt such that |m| < Mt 
and \v\ < Mt in [0,r] x [0,oo). We cannot confirm the uniqueness of the solution to (j2.ip 
without the assumption of boundedness since v is defined in an unbounded domain. For our 
problem (jl.ip . the solution is always bounded, see Theorem 12.51 below. 

Remark 2.3. It follows from the uniqueness of the solution to ()2.ip and a standard compactness 
argument that the unique solution (u, v, h) depends continuously on the parameters appearing in 
()2.ip . This fact will be used in the sections below. 

Theorem 2.4. Under the assumptions of Theorenn \2.1\ if we assume further that there exists a 
constant L > such that f{u, v) < L[u + v) and g{u, v) < L{u + v) for u,v > 0, then the unique 
solution obtained in Theorem \2.1\ can he extended uniquely to all t > 0. 

Proof: Let [0,Tmax) be the maximal time interval in which the solution exists. By Theorem 
12. H Tmax > 0. It remains to show that T^ax = 00. 

Suppose for contradiction that T^ax < 00. Fix M* G {Tmax, 00). Let {U{t),V{t)) be the 
solution to the following ODE system: 

.2l3^ ( Ut = LiU + V), Vt = L{U + V), t > 0, 

1 f^(0) = lko||c([o,M)' ^(0) = ||wo||l->([o,oo))- 

It is easy to see that 

<U + V < {\\uo\\c[o,ho] + lbo||L-'([o,+oo)))e^^* < (PollcicM + Po||L°°([o,+oo)))e^'^^* 
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for t G [0, Tmax)- Recalling the assumption that f{u,v) < L{u + v) and g{u,v) < L{u + v) and 
comparing {u,v) with {U,V) yield that for t G [0,Tmax) and r G [0, 

< u{t, r) + vit, r) < U{t) + V{t) < (||no|b[o,M + ll^ollL-([o,+oo)))e'^'^* := Cq. 

Next we claim that < h'{t) < Cj for all t G (0,Tjnax) and some C^ independent of Tmax- In 
fact, by the strong maximum principle and Hopf boundary lemma h'{t) is always positive as long 
as the solution exists. To derive an upper bound of h'{t), we define 

n = ^^M := {{t, r) -.OKtK T^^, h{t) - M'^ <r < h{t)} 

and construct an auxiliary function 

u{t,r) := Ce[2M{h{t) - r) - M^{h{t) - rf]. 

We will choose M so that u{t, r) > u{t, r) holds over J7. 
Direct calculations show that, for {t, r) G $7, 

ut = 2C6Mh'{t){l - M{h{t) - r)) > 0, 
-Ur = 2C6M[1 - M{h{t) - r)] > 0, 
_ _ 7V-1_ 2 

— Au = —Urr Ur > 2C%M , 

r 

f{u,v) < L{u + v) < LCq. 

It follows that 

Ut - diAu > 2di CeM^ > f{u, v) in ft 
if > 25^. On the other hand, 

u{t, h{t) - M-^) = Cq> u{t, h(t) - M"^), u{t, h{t)) = = u(t, h{t)). 

To use the maximum principle over $7, we only have to find some M independent of Tmax such 
that UQ{r) < u{0,r) for r G [/io — M^^, ho]. We calculate 

Ur{0,r) = -2CqM[1 - M{ho - r)] < -CqM for r e [ho - (2M)-\/io]. 

Therefore upon choosing 



M := max ■ 

we will have 



L M\uo\\c^([o,ho]) 
2di ' 3C6 



Ur{0,r) < n[,(r) for r e [ho - {2M)-^,ho]. 
Since u{0, ho) = Uo{ho) = 0, the above inequality implies 

u{0,r) > uo{r) for r e [ho - {2M)-^,ho]. 

Moreover, for r e [ho - M~^,ho - (2M)"^], we have 

3 3 
u{0,r) > -Ce, uo(r) < hoWcmoM)^"^ - 

Therefore no(r) < u(0,r) for r G [/lo - /io]- 

Applying the maximum principle to u — n over Q. gives that u{t,r) < u{t,r) for {t,r) G 
which implies that 

Ur{t,h{t)) > Ur{t,h{t)) = -2MCq, h'{t) = -nur{t,h{t)) < Ct := 2MCen 

for t G [0,rniax)- Moreover, 

ho < h{t) <ho + C-jt <ho + C7M*. 
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We now fix 5q £ {0,Tmax)- By standard parabolic regularity, we can find Cg > depending 
only on M* , L,Cq and C7 such that Ollci+^Clo./iW]) < C's and Ollci+^dCoo)) < C's for 
t S [^Oj^max)- It then follows from the proof of Theorem 12.11 that there exists a r > depending 
only on M*,L, Cq, C7 and Cg such that the solution of problem ()2.ip with initial time T^ax 
can be extended uniquely to the time T^ax — t/2 + r. This contradicts the maximality of T^ax- 

□ 

We have the following estimates. 

Theorem 2.5. Problem (II. ip admits a unique and uniformly bounded solution {u,v,h). That is, 
the solution is defined for all t > and there exist constants Mi and M2 such that 

< u{t, r) < Ml for t G (0, +00), < r < h{t), 

< v{t, x) < M2 for t £ (0, +00), < r < +00. 
Moreover, there exist a constant M3 such that 

< h'{t) < M3 for t e (0, +00). 
Further more, (jl.ip does not have any unbounded solution. 

Proof: By Theorem 12.41 (jl.ip has a unique bounded solution defined for all t > 0. It follows 
from the comparison principle that u{t,r) < u{t) for t € (0, 00) and r S [0,/i(t)], where 

ai iLt/ , , ai 

u{t) := 7— e*"! I e''i - 1 + 



hi V 6i||uo 
which is the solution of the problem 



(2.14) J -^ = u{a'i-biu), t>0, 



u(0) = ||no| 



Thus we have 



u{t,r) < Ml := supM(t). 

t>0 



Since v(t,r) satisfies 

( vt — d2^v < v{a2 — C2v), t > 0, < r < 00, 
I v{0,r) = vo{r) > 0, < r < 00. 

we have v{t,r) < max{||uo||L°o(o,+oo), ^} - M2. 

Using the strong maximum principle to the equation of u we immediately obtain 

u{t, r) > 0, Ur{t, h{t)) < for t > 0, < r < h{t). 

Hence h'[t) > for t € (0, 00). Similarly we have v{t, r) > for t > 0, < r < 00. 

It remains to show that h'{t) < M3 for t G (0, +00) and some M3. The proof is similar as that 
of Theorem 12.41 with Cq replaced by Mi and M3 = Cj = 2MMifi, we omit the details. 

We next show that any solution of (jl.ip is bounded, namely, there exists M > such that 
u, V < M in the range they are defined, whenever (u, v, h) is a solution to (jl.ip defined in some 
maximal interval t £ (0,T). Indeed, let U{x) be the unique boundary blow-up solution of 

-diAU = U{ai - biU) in Bi(0) := {r < 1}, U = 00 on a5i(0), 

and denote u{t,x) = |x|); then it is easily checked by using the comparison principle that 
u{t,xo + x) < ||uo||oo + U{x) for X G Bi{xo) and t > 0. It follows that u < ||uo||oo + U{0) in 
the range that u is defined. Similarly we can show v < ||fo||oo + ^(0), where V{x) is the unique 
boundary blow-up solution of 

-da Ay = V{a2 - C2V) in Bi{0), V = 00 on dBi{0). 

(The existence and uniqueness of U and V is well known; see, for example, [6].) □ 
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In what follows, we discuss the comparison principle for (jl.ip . For a given pair of functions 
u := {u,v) and u := {u,v), we denote 

[u,u] = {u ■= {u,v) e [C([0,r] X [0,oo))]2 : {u,v) < {u,v) < {u,v)}, 

where by (ui , wi ) < (n2, ^2), we meaniii < U2 and < ^2- A function pair (/, 51) = {f{u,v),g{u,v)) 
is said to be quasimonotone nonincreasing if for fixed n, / is nonincreasing in v, and for fixed v, 
g is nonincreasing in u; this is satisfied by / = u{ai — biu — civ) and g = f (02 — — C2v) in 
(fTTT) for li, i; > 0. 

Lemma 2.6 (The Comparison Principle). Let {f,g) be quasimonotone nonincreasing and Lips- 
chitz continuous in [u,u], with f{0,v) = g{u,0) = 0. Assume that T G (0, 00), h,h £ C^{[0,T]), 
ue C(D^)r\C^'^{D^) withD^ := {(t,r) gR"^ :tG {0,T],r e {0,h{t))}, u £ C(D^)nC^'^{D^*) 
with := {{t,r) G : t G (0,r],r G [0,7i(t))}, v,v£ (L°° n C)([0, T] x [0, 00)) n Ci'2((0, T] x 
[0,00)) and 



ut — diAu > f{u 








0<t<T, 


< r < h{t), 


Uf — diAu < f{u 


v), 






0<t<T, 


< r < hit), 


vt - d2Av > g{u, 


v), 






<t<T, 


< r < 00, 


Vt - d2Av < g{u, 


v), 






0<t<T, 


< r < 00, 


Ur{t,0) = V,.(t,0) 


= 


u{t,r) 


= 0, 


<t<T, 


h{t) < r < 00 


U^lt,0) =Vr{t,0) 


= 


u{t, r) 


= 0, 


< t < T, 


h{t) < r < 00 


h!{t) < -mrit^Kt)), 


h'it) > 


-fIUr{t,h{t)), 


< t < T, 




MO) <ho< h{0) 












u{0, r) < uo{r) < 


u(0, 


r), 




0<r <ho 




v{0, r) < VQ{r) < 


U(0, 


r), 




< r < 00 





(2.15) 



Let {u,v,h) be the unique bounded solution of (j2.ip . Then 

h{t) <h{t) in (0,r], u{t,r) <u{t,r), v{t,r) >v{t,r) for {t,r) G (0,r] x [0,oo), 
h{t) >h{t) in (0,r], u{t,r) >u{t,r), v{t,r) <v{t,r) for {t,r) G (0,r] x [0,oo). 

Proof: We only prove u < u, v > v and h < h; the result involving {u,v,h) can be proved 
in a similar way. Let M be an upper bound of v and v in [0,T] x [0, +00), w = M — v and 
w = M — V, then (u, W, h) satisfies 

' Ut - diAu > f{u, M -w), < t < T, < r < h{t), 

wt - d2Aw > -g{u, M -w), <t <T, < r < 00, 

Urit, 0) = Writ, 0) = 0, u{t, r) = 0, < t < T, h{t) < r < 00, 

-fMr{t, h{t)) < h'{t), <t<T, 

ho < h{0), uo{r) < u{0, r), < r < ho, 

M - vo{r) <w{0,r), < r < 00. 

First assume that ho < h(0). We claim that h{t) < h{t) for all t G (0, T]. If our claim does 
not hold, then we can find a first t* < T such that h{t) < h{t) for t G (0,t*) and h{t*) = h{t*). 
It follows that 

(2.16) h'{t*) > h'{t*). 

We now show that {u,w) < {u,w) in [0,t*] x [0,oo). Letting U = {u — u)e~^^ and W = 
{w — w)e~^^, we obtain 

' Ut - diAU > {~K + bii)U + bi2W, 0<t<t*, 0<r < h{t), 

Wt - d2AW > b2iU + {-K + b22)W, 0<t<t*, 0<r<oo, 

Ur{t,0) = Wr{t,0) = 0, U{t,r) = 0, 0<t<t*, h{t) < r < 00, 

C/(0, r) > 0, 1^(0, r) > 0, < r < 00, 



(2.17) 
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where K is sufficiently large such that K > 1 + \bii\ + 612 + I622I + &21 in [0, i*] x [0,oo), 
bii{t,r),b22{t,r) are bounded, b2i{t,r),bi2{t,r) are bounded and nonnegative since (/, 5) in (12.15P 
is quasimonotone nonincreasing and Lipschitz continuous. 

Since the first inequality of (|2.17p holds only in part of [0, 00), we cannot use the maximum 
principle directly. We first prove that for any / > h{t*), 

rr, N M(r^ + 2dNt) , M(r^ + 2dNt) 
U{t,r) > ^ ^ and W{t,r) > ^ ^ 

in [0, t*] X [0, /], where d = max((ii, ^2). (Recall that M is an upper bound of v and v-) 

We observe that due to the inequalities satisfied by u, we can apply the maximum principle to 
u over the region {(t, r) : < r < h{t), < t < T} to conclude that n > 0. Set 

I7(t, r) = u+ and Wit, x) = W + 

then due to our choice of K, (U, W) satisfies 

Ut - diMJ_^ i-Kj- bii)U + bi2W^ 0<t<t*, 0<r < h{t), 

Wt-d2AW>b2iU+i-K + b22)W, 0<t<t*, 0<r<l, 

Ur{t,0) = Wr{t,0) = 0, U{t,r) > > 0, 0<t<t*, h{t) <r<l, 

W{t,l) = W(t,l) + ^'^(''+^'^^*) > 0, 0<t<t*, 
I U{0, r) > 0, W{0, r) > 0, < r < /. 

We now prove that min{min[o^f*]x[o,«] ™iii[o,t*]x[o,i] ^} ■= t > 0. In fact, if r < 0, then there 
exists (io)^o) S with < to < t* and < tq < /i(to) such that C/(to,To) = r < 0, or there 
exists (ii,ri) G with < ti < t* and < ri < / such that W{ti,ri) = r < 0. For the former 
case, (Ut - diAU){to,ro) < 0, but 

[i-K + bu)U + bi2W]ito,ro) > {-K + \bu\)T + bi2T>-T>0. 

For the latter case, (Wt - d2AW){ti,ri) < 0, but 

[b2iU + i-K + b22)W]{ti,ri) > {-K + \b22\)T + b2iT > -r > 0. 

Both are impossible. Therefore r > 0, that is U > and > in [0,t*] x [0,/], which implies 
that 

, M{r'^ + 2dNt) , M{r^ + 2dNt) 
U{t,r)> ^-^2 W{t,r)> ^-^2 

iov <t <t*, <r <l. Taking / 00 yields that U{t, r) > and W{t, r) > in [0, t*] x [0, 00), 
therefore u <u and w <w in [0, t*] x [0, 00). 

We now compare u and u over the bounded region 

:= {{t,r) G : < t < t*,0 < r < h{t)}. 

Since Z(t, r) := u(t, r) — u{t, r) satisfies 

Zt - diAZ > bnZ + buiw - w) > buZ, 0<t<t*, < r < h{t), 

the strong maximum principle and the Hopf boundary lemma yield Z(t,r) > in fi^*, and 
Zr{t*,h{t*)) < 0. We then deduce that h'{t*) < h'{t*). But this contradicts (l2T6]l . This proves 
our claim that h{t) < h{t) for all t G (0,T]. We may now apply the above procedure over 
[0,T] X [0, oo) to conclude that u<u and w <w (i.e. v>y) in [0,T] x [0,oo). Moreover, u <u 
for t G (0,r] and r G [0,h{t)). 

If ho = /i(0), we use approximation. For small e > 0, let {ue,Ve, he) denote the unique solution 
of (jl.ip with /iQ replaced by /io(l — e). Since the unique solution of (jl.ip depends continuously on 
the parameters in (jl.ip . as e — 0, (ne,Ue,/ie) converges to {u,v,h), the unique solution of (II. ip . 
The desired result then follows by letting e — )• in the inequalities <u, > v and < h. □ 
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Remark 2.7. The conclusions in Lemma 12.61 remain valid if the condition 

Zr{t,0) = for z G {u,u,v,v} 

is replaced by 

Zr{t, 0) > for z e {u,v}, Zr{t, 0) < for z € {u,v}- 

To see this, we only need to observe that all the arguments in the proof carry over to the new 
case, except that we now have to avoid tq = or ri = 0, since the functions U{t,\x\) and 
W{t, \x\) are not in x at x = 0. However, the above conditions guarantee that Ur{t,0) < 
and Wr{t,0) < 0. Therefore (t,0) cannot be a minimum point of these functions. This implies 
that ro > and ri > 0. 

We next fix UQ,VQ,di,ai,bi,Ci and examine the dependence of the solution on fi, and we write 
{u^,v^, h^) to emphasize this dependence. As a consequence of Lemma 2.6, we have the following 
result. 

Corollary 2.8. For fixed Uo,VQ,di,ai, hi andci. If < ^2- Then u^^ {t,r) < u^^{t,r), v^^{t,r) > 
v^^{t,r) fort G (0,oo), r G [0,h^'^{t)) and h^'^it) < h>^^{t) in (0,oo). 

3. Invasion of an inferior competitor 

In this section, we examine the case that u is an inferior competitor, namely 

/o IN ai . J^i ci 

(3.1) — <mm<— , — 

The following theorem shows that the inferior invader cannot establish itself and the native 
species always survives the invasion. 



Theorem 3.1. // (13. 1|) holds and vq ^ 0, then 

lim (u(t, r), t;(t, r)) = (0, — ) uniformly in any compact subset o/[0,oo). 
i^+oo \ C2 J 

Proof: First we recall that the comparison principle gives u{t,r) < u*{t) for t > and 
r G [0, h{t)], where 

M (t) = — e''i e*"! + — — - — 
oi ^ oill-u-olloo 

is the solution of the problem 

{u*y = u*{ai-biu*), t>0, 
u*{0) = \\uo\\oo- 



(3.2) 

Since limt^oo u* (t) = we deduce 



lim sup n(t,r) < ^ uniformly for r G [0, oo). 

Similarly, we have 

(3.3) limsupu(t, r) < — uniformly for r G [0,oo). 

t^+oo C2 

Therefore for ei = (f|-2i)/2, there exists ti > such that u(t, r) < f^ + ei fort > ti,r G [0,oo). 
Then v satisfies 

Vt — d2Av > v{b2£i — C2v), t > ti, < r < oo, 

(3.4) { Vr{t,0)=0, t>ti 

v{ti,r) > 0, < r < oo. 
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Let f^, be the unique solution to 



{v*)t — d2/S.v^, = v^{b2ei — C2f*), t > ti, < r < oo, 
(v*V(t,0) = 0, t>ti 
v*{ti,r) = v{ti,r), < r < oo. 

It is well known ([6j) that limt^oo v<t{t,r) = 62^1 /c2 uniformly in any bounded subset of [0,oo) 
Therefore for any L > 0, there exists ti > ti such that 



(3.5) 



v{t, r) > r) > 



2C2 



for t > iL, < r < L. 



Now (u, v) satisfies 



(3.6) 



ut — diAu = u{ai — biu — civ), t > ti, < r < h{t), 

vt — d2Av = v{a2 — b2U — C2v), t > t^, < r < 00, 

Ur{t,0) = Vr{t,0) = 0, t > ti, 

[ u{t,r) <^ + ei, v{t,r) t>tL, 0<r<L. 



Since u = for t > tL, r > h(t), no matter whether or not h(t) < L, we always have u < u and 
V > in [ii,oo) x [0,L], where {u,v) satisfies 



(3.7) 



Ut — diAu = u{ai — biu — civ), t > t^, < r < L, 
v_t — d2A.v_ = v{a2 — b2U — C2v), t > tL, < r < L, 

Ur{t,0) =Vr{t,0) = 0, 

. u{t,r) = '^ + ei, v{t,r) 



t > tL, 

t>tL, r = L,OTt = tL,G<r<L. 



The system ()3.7p is quasimonotone nonincreasing, which generates a monotone dynamical system 
with respect to the order 

{ui,vi) <p {u2,V2) if and only if ui < U2 and vi > V2, 

with the initial value + ei, an upper solution. It follows from the theory of monotone 

dynamical systems (see, e.g. Corollary 3.6) that limt_s>+oo u{t, r) = UL{r) and linit^+oo v_{t, r) 
u^(r) uniformly in [0,L], where {ul,IIl) satisfies 



(3.8) 



-diAuL = UL{ai - biUL - ciVl), <r < L, 
-d2AvL = VL{a2-b2UL - C2Vl), 0<r<L, 



dr dr 



^-^(o) = l^(o) = o 



[ UL{L) = ^ + e,, vl{L) = '-^, 



and is the maximal solution below (^ + ei, ^^) of the above problem under the order <p. 

Next we observe that if < Li < L2, then UL^{r) > ul2{^) aiid VL^^ir) < VL2i''') 
which can be derived by comparing the boundary conditions and initial conditions in ()3.7p for 
L = Li and L = L2. 

Let -L — )• 00, by classical elliptic regularity theory and a diagonal procedure, it follows that 
{uL{T),v_L{r)) converges uniformly on any compact subset of [0, 00) to iuoo-,v_^), which satisfies 

-diAHoo = Uooiai - biUoo - civ^), < r < 00, 

-d2Avoo = Vooia2 - huoo - C2V00), < r < 00, 
9f„ 



^(0) = ^(0) = 0, 



< r < 00. 
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Next we show that (r) = and Voo{r) = ^. To this end, let us consider the fohowing ODE 
system: 

{Zt = z{ai — biz — ciw), t > 0, 

Wt = w{a2 - h2Z - C2w), t>0, 
ziO) = ^^ + eu ^0) = ^. 

Since ai/a2 < min{ci/c2, 61/62}, it is well-known (e.g.[T2]) that {z,w) — ?■ (0, as t — >• 00. 
Therefore the solution (Z{t,r), W{t,r)) of the problem 

' Zt - diAZ = Z{ai - hiZ - ciW), t > 0, r > 0, 
Wt - d2AW = W{a2 - b2Z - C2W), t > 0, r > 0, 

Zr{t,0) = Wr{t,0) =0, t>0, 

. Z{0,r) = ^ + ei, W{0,r) = '-^, r>0 

satisfies {Z, W) — )■ (0, ^) as t — )• 00 uniformly in [0, 00). It follows from the comparison principle 
that Uooir) < Z(t,r) and Vodr) > W{t,r) for t > 0, which immediately gives that = and 



(3.10) 



V =22 



We thus have limsupj_^^oo ^(^1 < and liminf(_^+oo v{t, r) > ^ and uniform in [0, L], which 
together with ()3.3p implies that limj_^+oo uit, r) = and limt_j.+oo I'it, ^) ~ % uniformly in any 
bounded subset of [0,cxd). □ 

Remark 3.2. The above proof can be used to show that the unique solution (u(t, j;), f (t, x)) 
of the Cauchy problem of (jl.Sp . with m(0, x), v(0, x) nonnegative, bounded and t;(0, x) ^ (but 
not necessarily radially symmetric), also converges to R2 = (0,02/02) locally uniformly in as 
t — )• 00 when (13.11) holds. 



We note that Theorem 13.11 gives no information on the dynamical behavior of the spreading 
front, and the exact behavior of {u,v) over the entire spatial range < r < 00 is also unclear. 
Our next result provides such information, provided that the native species is already rather 
established at t = 0, in the sense that infr>o VQ{r) > 0. 

Let us observe that Theorem 12.51 implies h{t) is monotonic increasing and therefore there exists 
hcxi £ (0,-|-oo] such that limt^+00 h{t) = h^o- 

Theorem 3.3. Suppose that ()3.ip holds and vo{r) > 5 > for < r < 00. Then hoo < 00 and 
limt_^+oo ^i(i, ?^) = 0, limt^+00 v{t,r) = ^ uniformly in [0,oo). 

Proof: First it follows from the comparison principle that 

limsupn(t,r) < — !- uniformly for r G [0, 00). 

bi 

limsupt'(t, r) < — uniformly for r G [0, 00). 

C2 

Therefore for £1 = (||-|^)/2, there exists ti > such that ii(t, r) < f^ + ei fort > ti,r e [0,oo). 

Note that < u{t,r) < Mi := max{2i, ||mo||c([o,/io])} and < v{t,r) < M2 := max{^, ||wo||loo([o,+oo))j 
for t>0,0<r<oo. Therefore v satisfies 

Vt - d2Av > v{a2 - b2Mi - C2M2), t > 0, < r < 00, 
Vr(t,0)=0, t>0 

v{0,r) > 6, < r < 00, 

and hence v{t,r) > Je^-^^Afi-caMa)* for t > and < r < 00. 
Now let us consider the following problem: 

Zt = z{ai — biz — ciw), t > ti, 

(3.11) < Wt = w{a2 - b2Z - C2w), t > ti, 
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It follows from the comparison principle that u{t,r) < z{t) and v{t,r) > w{t) for t > ti,0 < r < 
oo. Under the assumption 01/02 < min{ci/c2, 61/62}; it is well-known that {z,w) — (0, fj) as 
t —7- 00. It follows that limt^_|_oo u{t, r) = uniformly for r G [0, 00) (we note that n(t, r) = for 
r > h{t)). Next we prove limt^oo v{t,r) = ^ uniformly for r G [0, 00). Since limj_!.oo ti(t, r) = 
uniformly for r G [0, 00), for any e > there exists T > such that < u{t, r) < e for t>T and 
r G [0,00). Recall that v{t,T) > Se^-^^^h-c2M2)T _ ^Ym^ ^^^^ 

vt — d2^v > v{a2 — 626 — C2v), t >T, < r < 00, 
Vr{t,0)=0, t>T, 
v{T, r) > 5e(-^2Afi-c2Af2)T^ < r < 00 

Consider the following problem 

Vt = v{a2 - b2e - C2v) for t > T, v{T) = Je(-^2Afi-c2M2)T^ 

It follows from the comparison principle that v{t) < v{t,r) for t > T and r G [0,oo). Since 
'^{'t) ~^ "'2-£b2 ^ _^ have "-^"^^'-^ < liminf(_>>oo f (t, r) uniformly for r G [0,oo). Since 

e > is arbitrary we have limt^c>ov{t,r) = ^ uniformly in [0,oo). 

Next we prove that h^o < 00. Since limt^^oo u(t,''') = and limt^^oo v(t,r) = ^ uniformly 

'-■2 

in [0, 00), we have limf^+oo(oi — biu{t, r) — cif (t, r)) = oi — 02C1/C2 < uniformly in [0, cxd), and 
there exists T* such that (oi — biu{t, r) — civ{t, r)) < for i > T*, < r < 00. 
Direct calculation gives 

J fh{t) fh{t) 

^ / r^~\{t,r)dr = / r^-\tit,r)dr + h^-\t)h'{t)uit,hit)) 
di Jo Jo 

rh(t) ph(t) 



I dir^^-'^Andr + / n(ai — 6in — civ)r^~^dr 
Jo Jo 
ri rh(t) 
-h^-\t)h'{t) + / u{ai-biu-civ)r^-^dr. 



Integrating from T* to t yields 

V-l /, N 1 / ,rT.* , "1 , V/^*\ "1 



(3.12) 



0< / r^-in(t,r)dr= / r'''^uiT\r)dr + ^h"" (T*) - ^h"" it) 
Jo Jo . . ^/^ 



i-t Ms) 

+ I I u{ai — biu — civ)r'^^^drds, t>T*. 
Jt* Jo 



Since (oi — biu{t, r) — civ{t, r)) < for t > T*, < r < 00, we deduce 

h^{t) < -f / r^-\iT*,r)dr + h^{T*), t > T*, 
"1 Jo 

which implies that /iqo < 00. □ 
Theorem 13.31 suggests that an inferior competitor can never penetrate deep into the habitat of 
a well established native species, and it dies out before its invading front reaches a certain finite 
limiting position. 

4. Invasion of a superior competitor 

This section is devoted to the case that u is a superior competitor, that is 

u -i\ ai J61 ci 

(4.1) — > max < — , — 

a2 [ 62 C2 

Let Ai(i?) be the principal eigenvalue of the operator —A in Bji subject to homogeneous 
Dirichlet boundary conditions. It is well-known that Ai(i?) is a strictly decreasing continuous 



(4.3) 
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function and 

lim Xi{R) = +00 and lim Xi{R) = 0. 

Therefore, there exists a unique R* such that 

(4.2) Xi{R*) = 1. 

It is easy to check that R* = 7r/2 if = 1. 

For the sake of convenience and completeness, we first recall the following spreading-vanishing 
dichotomy for the radially symmetric diffusive logistic problem 

ut — dAu = u(a — bu), t > 0, < r < h{t), 

Ur{t, 0) = 0, u{t, h{t)) =0, t > 0, 

h'{t) = -nurit,h{t)), t>0, 

^ h{0) = ho,uiO, r) = uo{r), < r < /iq. 

Proposition 4.1. Let (u(t, r), /i(t)) be the solution of the free boundary problem (|4.3p . Then the 
following alternative holds: 
Either 

(i) Spreading: hoo = +00 andlimt^Ji-oou{t,r) = | uniformly in any compact subset o/[0,oo); 
or 

(ii) Vanishing: hoo < R* \fi and limf^+oo \\u{'t,-)\\c{[ii,h(t)]) = ^■ 

Proposition 4.2. // /iq > R*\[\' then spreading always happens. If ho < R*\l~ji, then there 
exists n* > depending on uq such that vanishing happens when n < fj,* , and spreading happens 
when fj. > fi* . 

The proofs of Propositions 14. l l and 14.21 can be found in [5j for the one dimensional case and [3] 
for higher space dimensions. 

We will show that when ()4.ip holds, a similar spreading-vanishing dichotomy holds for p.ip . 
More precisely, we have the following results. 

Theorem 4.3. Suppose that (j4.ip holds and (u, v, h) is the unique solution of (II. ip with vq ^ 0. 
Then the following alternative holds: 
Either 

(i) Spreading of u: h^o = +00 and Ymit^j^oa {u{t,r),v{t,r)^ = 0^ uniformly in any 
compact subset of [0,oo); 

or 

(ii) Vanishing of u: h^o < R* \/ ai-a2c{ J^ ^'^^ limt^+oo hit, ■)\\c{[oMm = 0, limt^oo v{t,r) = 
^ uniformly in any compact subset of [0, 00). 



Theorem 4.4. In Theorem \J^.3\ if h^ > R* y ai~a2Ci/c2 ' then spreading of u always happens. If 
ho < R* \J '; ^a2Ci/c2 ' then there exists fi* G [0, 00), depending on (uo,vo), such that spreading of 
u happens exactly when fi > fi* . Moreover, ^u* > when ho < R* ^J^. 

We prove these results by several lemmas. In the rest of this section, we always assume that 
(j4.ip holds, and {u,v,h) is the unique solution of (jl.ip . with vq ^ 0. 

Before starting the proofs, let us note that by the symmetric positions of u and v in (jl.Sp . 
we may use Remark 13.21 to conclude that when (14. ip holds, the unique solution (u, v) of the 
Cauchy problem of (jl.Sp with u{0, x),v{0, x) bounded, nonnegative and not identically zero. 
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satisfies limj_j.oo {u{t,x),v{t,x)') = Ri = y^, Oj locally uniformly in M^. The behavior of the 

free boundary problem (jl.ip described by Theorems 14.31 and 14.4^ however, is more complicated, 
where vanishing for u is possible. 

Lemma 4.5. If hoc = oo, then limt^^oo u{t,r) = |^ and limt^^oo v{t,r) = uniformly in any 
bounded subset of [0, oo). 

Proof: First we recall that the comparison principle gives 

limsupu(t,r) < -i uniformly for r G [0, oo), 



02 



limsupf(t,r) < — uniformly for r S [0, oo). 

t-^+oo C2 

Therefore for 82 = - ^)/2, there exists ^2 > such that v{t, r) < ^ + £2 for t > t2, r e [0, 00). 
Therefore u satisfies 

ut — di/S.u > u{cie2 — biu), t > t2, < r < h{t), 
Ur{t, 0) = 0, u{t, h{t)) =0, t> t2, 

h'{t) = -flUr{t, h{t)), t > t2, 

^ u{t2,r) > 0, < r < h{t), 



(4.4) 



Next, we prove that for any I > R*^/ -^r^, there exists ti > ^2 such that 



(4.5) 



u{t, r) > 



Cl£2 

2bi 



for t > t/, < r < I. 



Indeed, since h{t) — )• 00, there exists is > t2 such that h{t^) > I. It follows from the comparison 
principle ([!]) that u{t,r) > Ui{t,r) for t > t^, r £ [0,h{t)] and h{t) > h{t) in (t3,oo), where 
{uijh.it)) is the solution of the following free boundary problem discussed in |4]: 

{ui)t - diAu; = U;(cie2 - biui), t > ts, < r < h{t), 
{ui)r{t,0)=uiit,h{t)) = 0, i>*3, 

U{t) = -mr{t,h{t))i t>ts, 

[ lkit3,r) = u{t3,r), hits) = hits), < r < hits). 



(4.6) 



Since cie2 > cii(-^)^, it follows from Proposition 14.21 that hit) — )• 00 and Uiit,', 
uniformly in [0,/], which implies that ()4.5p holds for some ti. 
Now we know that (u, v) satisfies 



^ast 



00 



(4.7) 



Ut — diAu = ti(ai — biu — civ), t > ti, < r < I, 
vt — d2^v = f (02 — b2U — C2v), t > ti, < r < I, 

Urit,0) = Vrit,0) = 0, t>tl, 

[ uit,r) > q§^, vit,r) < f +£2, t>ti, 0<r<l. 

As in the proof of Theorem 13.11 it follows from the theory of monotone dynamical systems that 
liminfj_).+oo ""(t, r) > u^ir) and limsupt_^_^^^vit,r) < Viir) in [0,/], where iui,vi) satisfies 

—diAui = Uiiai — biUi — civi), < r < I, 
-d2Avi = viia2 - b2'Ua - C2V1), < r < /, 



du 
dr 

[ niil) 



= ^(o) = o, 
^> Ml) - 



C2 



+ £2- 



Letting I — )• 00, we similarly have iui,vi) — )• (n^ 



which satisfies 



-diAu^ = Moo(ai - huoc 
-d2AVoo = Vooia2 - huoo ■ 



- CiVoo), 
C2V00), 



^(0) = 0, 



< r < 00, 
< r < 00, 

< r < 00. 
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Due to (j4.ip . as before, by the global dynamical behavior of the associated ODE system 
([T2]). we deduce that u^{r) = ^ and Voo{f) = 0. We thus obtain liminft^+oo ''") > ^ 
and Iimsup^^_|_oo r) < uniform in [0,/], which implies that limt-^^oou{t,r) = ^ and 
limj_>.+oo v{t, r) = uniformly in any bounded subset of [0, oo). □ 

Lemma 4.6. Ifhoo < oo, i/ien limt^+oo \\uit, ■)\\c([0Mt)]) = ^ andlimt^+oov{t,r) = ^ uniformly 
in any bounded subset of [0,oo). 



Proof: Define 

' ~ Kt) 

Direct calculations yield 



w(t, s) = u(t,r), z(t, s) = v{t,r). 



Hence w{s,t) satisfies 



(4.8) 



wt - dij^AsW - ^sws = w{ai - biw - ciz), < s < ho, t > 0, 
Wsit,0) = w{t,ho) = 0, t>0, 
w{0, s) = no(s) > 0, < s < /iQ. 



\ai — biw — ciz\\loo < oi + 6iMi + C2M2, 



<M3. 



This is an initial boundary value problem over a fixed ball {s < /iq}- Since ho < h(t) < h^o < 00, 
the differential operator is uniformly parabolic. By Theorem 2.5, 

m 
hit)- 

Therefore we can apply standard theory and then the Sobolev imbedding theorem ([10]) to 
obtain, for any a S (0, 1), 

ll^llc(i+°)/2.i+"([0,oo)x[0,ho]) - 

where C4 is a constant depending on a, /iq. Mi, M2, M3 and ||mo||ci+"[o,/io]- Similarly we may use 
interior estimates to the equation of z to obtain 

||2:|lc(i+a)/2,i+«([o,oo)x[0,ho]) - ^5, 

where C5 is a constant depending on a, /io, Mi, M2, M3 and ||i'o||ci+"[o,ho+i]- 

It follows that there exists a constant C depending on a, /iq, (wq, wq) and h^ such that 

(4.9) ||n||p(H-Q)/2,l + a(G) + ||t'||c(l + c)/2,l + <:«((J) + 1 1 ^ 1 1 C1 + q/2 ( [0,oo)) — , 

where G := {(t,r) : t > 0,r G [0,h{t)]}. 

Arguing indirectly, we assume that lim sup(_j._|_oo ||'u(t, •)llc([o,/i(i)]) = 5 > 0. Then there exists 
a sequence {tk-,rfS) with < < cxd, < r^. < h{tk) such that u{tk-,rk) > 5/2 for all /c G N, 
and tfc — )• 00 as A; — )• 00. Since u{t,h{t)) = and since (j4.9p infers that \ur{t,h{t))\ is uniformly 
bounded for t S [0, 00), there exists o" > such that < h{tk) — a for all k >1. Therefore a 
subsequence of {r^} converges to tq S [0, hoc — cr]. Without loss of generality, we assume — )• vq 
as — >• 00. 

Define 

Uk{t, r) = u{tk + t, r) and Vk{t, r) = v{tk + t, r) for (t, r) G Gfc 

with 

Gfc := {{t, r):te {-tk, 00), r G [0, /i(tfc + t)]}. 
It follows from (14. 9p that {{uk,Vk)} has a subsequence {(^ifci, ^^fcj} such that 
\\iuk,,VkJ - (n,t5)||ci.2(GfeJxCi.2(GfcJ ^ as « ^ 00, 
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and {u,v) satisfies 

ut — diAu = u{ai — biu — civ), t G (—00, 00), < r < /too, 
vt — d2Av = v{a2 — — C2v), t G (—00, 00), < r < h^, 

with 

hoo) = for t G (—00, 00). 

Since n(0, rg) > (5/2, the maximum principle infers that -S > in (—00,00) x [0, /Iqo)- Thus we 
can apply the Hopf boundary lemma to conclude that gq := {2^.(0, hoo) < 0. It follows that 

Ur{tki,h{tkJ) = drUk^{0,h{tkJ) < ao/2 < 

for all large i, and hence h'{tk^) > > for all large i. 

On the other hand, since ||^||ci+«/2([o 00)) — ^'(*) > ^ ^(^) — ^oo' ^® necessarily have 
^ as t — >■ 00. This contradiction shows that we must have 

^Um^||u(t,-)||c([o,Mt)]) =0- 

We may now use a simple comparison argument to deduce that limt^-^oov{t,r) = ^ uniformly 
in any bounded subset of [0, 00). □ 

Lemma 4.7. If h^o < oo, then /too < / • Hence hn > R* \ , implies 

J °" > °" — Y ai— 02Ci/c2 " — Y ai— a2Ci/c2 ^ 

hoo = OO- 

Proof: Assume for contradiction that R*\/ ^^—7— < hoo < 00. Then there exists Ti > 

Y ai— a2Ci/C2 



such that I := h{Ti) > R* y ai-a2ci/c2 - ^® choose e sufficiently small such that h(Ti) > 



R* 



di 



ai—a2Ci/c2—ci£ ' 

Recall that 

limsupf(t, r) < — uniformly for r G [0, 00). 

t-^+00 C2 

Therefore for the above chosen e > 0, there exists T2 > Ti such that v{t, r) < ^ + £ for 
t>T2,r & [0,00). Hence {u,h) satisfies 

Ut — diAu > u{ai — Ci{a2/c2 + £) — biu), t > T2, < r < h{t), 

Urit,0) = 0, u{t,h{t)) =0, t>T2, 

h'{t) = -i^Ur{t,h{t)), t>T2, 

^ u{T2,r) > 0, 0< r < h{T2). 

This indicates that {u, h) is an upper solution to the problem 



(4.10) 



(4.11) 



Ut — diAn = u{ai — 01(02/02 + e) — biu), t > T2, <r < h{t), 

u^{t,0) = 0, u{t,h{t)) = 0, t>T2, 

Ilit) = -mrit,h{t)), t>T2, 

u{T2,r) = u{T2,r), h{T2) = h{T2), < r < h{T2). 



Hence h(t) > h{t) for t > T2. Since h{T2) > I > R* ■^J^^^^^j^^^, it follows from Proposition 

4.2 that h(t) 00 and therefore hoo = 00. This contradiction proves that hoo < ^* a^ a2ci/c2 ' 

□ 

Lemma 4.8. If h^ < R* '^^ a2ci/c2 ' ^^^'^ there exists fJ, > depending on uq and vq such that 
hoo = 00 when ijl> ji. 
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Proof: Since 



limsupT;(t, r) < — uniformly for r G [0, oo) 

t-s>+oo C2 



(4.12) 



for e = — fj)/2 > 0, there exists T, which is independent of such that v{t, r) < ^ + e for 
t>T,r £ [0,oo). Thus {u,h) satisfies 

ut — diAu > u{cie — biu), t > T, < r < h{t), 

Urit, 0) = 0, u{t, h{t)) =0, t>T, 

h'{t) = -nurit,h{t)), t>T, 

u{T,r) > 0, < r < h{T). 

Note that u{T,r) depends on /i, but u(T,r) > z(T,r) in [0, oo) x [0, /iq], where z{t,r) and w{t,r) 
are determined by 

Zf — diAz = z{ai — biz — ciw), t > 0, < r < ho, 

wt — d2Aw = w{a2 — b2Z — C2w), t > 0, < r < ho, 

Zrit,0) = Writ,0) =0, t>0, 

z{t,ho) = 0, w{t,ho) = max{f|, ||fo||L°°}, t > 0, 
z(0,r) = no(r), w{Q,r) = m.a.yi{^,\\vo\\L^}, < r < ho- 

Clearly, z{T, r) is independent of /U. Now it is easy to see that (n, h) is an upper solution to the 
problem 

Ut - diAu = u{cie - biu), t >T, <r < h{t), 
Ur{t,0) = 0, u{t,h{t)) = 0, t>T, 
h[{t) = -^iu^{t,h{t)), t>T, 
u{T, r) = z{T, r), h{T) = ho, < r < ho- 

By Lemma 2.8 of [1], there exists fi > such that h^ = +oo for fi > fi. One can actually argue 
as in Lemma 3.7 of [5] to show that if 



(4.13) 



^ di((R\ ^^^)^-C^ 

^ J 1 1 II ^H V "1 -'^2^1/^2'^ ° 

/i > /i := max < 1, \\z{^ , r) 



cie" ' "'~J N{f^°r^-^z{T,r)dr) ' 
then h^ = +oo. Therefore hoc = oo for /i > /u. □ 

Lemma 4.9. There exists /i* > depending on uo and vq such that /iqo = +oo if fi > fi* and 

hoo < oo if < n < fj-* ■ 

Proof: Define S:=|/i>0:/ioo> ^f^^^^ETJ^} ■ follows from Lemmas HZ] and SS] that 
fi* := inf S G [0,oo). By Lemma 14.71 and the monotonicity of h^o with respect to /u (Corollary 
12. Sp . we find that hoo = +oo when /i > and hoo < +oo when < fi < n*. 



We claim that if ^* > 0, then /u* S. Otherwise hoo > ^* y ai-~a2Ci/c2 = Hence we 

can find T > such that h(T) > ^* y ai-a2Ci/c2 ' stress the dependence of the solution (n, v, h) 



of (II. ip on fi, we now write (uu, Vn, hn) instead of (u, v, h). So we have hu* (T) > R* \/ _ \ i . 
By the continuous dependence of {Ufj_,v^, h^) on ji, we can find e > small so that h^(T) > 



y ai-a2ci/c2 ^ ^ [/"* ~ + f]- It fohows that for aU such fi, 

lim /i^(t) > h^{T) > R*J . 

t^oo V "1 ~ 0,2Ci/C2 

This implies that [^* — e, ^* + e] C S, and inf S < n* — e, contradicting the definition of /x*. This 
proves our claim that /i* G S. □ 
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Lemma [4 . 71 implies that the threshold constant ^u* = if /iq > R* 



ai—a2Ci/c2 



- — -, — . The next lemma 



shows that /x* > if /iq < R*J^- 



Then there exists > depending on uq such that 



(4.14) 



0, t>0, 
t > 0, 

0<r <ho. 



Lemma 4.10. Suppose ho < R* 
hoc < +00 if n <]!■ 

Proof: Clearly (u, h) satisfies 

ut — diAu < u{ai — biu), t > 0, < r < h{t), 
Ur{t,0) = 0, u{t,h{t)) 
h'{t) = -nur(t,h(t)), 
u{0,r) = Uo{r), 

This indicates that (n, h) is a lower solution to the problem 

Ut — diAu = u{ai — biu), t > 0, < r < h{t), 

Ur{t, 0) = 0, u{t, h{t)) = 0, t > 0, 

h'{t) = -ijMr{t,h{t)), t>0, 

{ u{0,r) = uo{r),h{0) = ho, 0<r<h{0). 

Since Hq < R*\/^, it follows from Proposition 14.21 that there exists 77 > depending on uq such 
that hoo < +00 if < 77, therefore /loo < +00 if n <]I. □ 
Theorems 14.31 and 14.41 now follow directly from the conclusions proved in the above lemmas. 

5. ESTIMATES OF SPREADING SPEED 

In this section we give some rough estimates on the spreading speed of h{t) for the case that 
spreading of u happens. We always assume that (j4.ip holds. 

We first recall Proposition 3.1 of jl], whose complete proof is given in [1]. 



(4.15) 



Proposition 5.1. For any given constants a>0,b>0, d>0 and k € [0, 2v ad), the problem 
(5.1) - dU" + kU' = aU - bU^ in (0, 00), U{0) = 

admits a unique positive solution U = Uk = Ua,b,d,k, and it satisfies U{r) ^ as r ^ +00. 
Moreover, Ul{r) > for r > 0, f/^^(0) > U'f.^{0), UkA^) > Uk^ir) for r > and ki < ^2, and 
for each fi > 0, there exists a unique k^ = ko{fi,a,b,d) G (0, 2\/ad) such that ^C/^^(0) = /cq- 
Furthermore, 



(5.2) 



lim -X = 2, lim 4^-^ = 1/^/3. 

ad ^^0\/adan 



Making use of the function kQ{fj,,a,b,d), we have the following estimates for the spreading 
speed of h{t). 

Theorem 5.2. Assume that (j4.ip holds. If h^o = +00, then 

A;o(/x, ai — a2Ci/c2, bi,di) < liminf — ^ < limsup — ^ < ko{fj,, ai, bi,di). 

t t->.+oo t 



Proof: Since 

Ut — diAu = u{ai — biu — civ) < u{ai — biu), t > 0,0 < r < h{t), 

Ur{t, 0) = 0, u{t, h{t)) = 0, t>0, 

h'{t) = -nur{t,h{t)), t>0, 

u{0,r) = Uo{r), < r < Hq, 
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the pair (n, h) is a lower solution to the problem 

ut — di An = u{ai — biu), t > 0,0 < r < h{t), 

Ur{t,0) =0, u{t,h{t)) = 0, t>0, 

h'{t) = -fMr{t,h{t)), t > 0, 

, u(0, r) = no(r), < r < ho, 

It follows that h{t) > h{t) — > oo as t — ^ oo. By [3], 



We thus have 



Next we prove that 



lim = ko{fJ-,ai,bi,di). 
r— ^oo t 

hit) , , , , , 

limsup < kQ[fi,ai,bi,dij. 

t—^+oo t 

koifx, ai — a2Ci/c2, foi, di) < liminf ■ ^ 



Since lim supf_j.oo v < ^ uniformly for r £ [0, oo) and h^o = oo, for any < e < Eq := — 

there exists T^, such that v(t,r) < ^ + e for t > T^, r G [0,oo) and h(Te) > R\/ „ .fS,. 

Hence (n, /i) satisfies 

tii — diAu > u{ai — Ci(a2/c2 + e) — ^'iit), t > T;,, < r < 
tir(t,0) = 0, = 0, i>r£, 

= -^Ur{t, h{t)), t > Ts, 

u(Te,r)>0, 0<r<h{Te). 

This implies that [u, h) is an upper solution to the problem 

u-i — diAu = u{ai — Ci(a2/c2 + e) — biu), t > T^, < r < h{t), 
Urit,0) = 0, u{t,h{t)) = 0, t>Ts, 

h^it) = -mr{t,h{t)), t>Ts, 
u{Te,r)=u{Te,r), 0<r<h{T,). 



(5.3) 



(5.4) 



By Proposition 14. H h^ = oo since h{T^) > R* \/ _ i , \ ■ Moreover, from [3] we have 



ai— ci(a2/c2+£) ' 

limj_>oo = ko{fj.,ai — 01(02/02 + e),bi,di), which implies that 

kQ{fj,,ai - ci (02/02 +e),bi,di) < liminft_>+oo ^ for any < e < Eq. 

Letting e — )• and using the continuity of kQ with respect to its arguments, we immediately 
obtain the desired result. □ 
Next we obtain an upper bound for the spreading speed that is independent of fi, namely, we 
show that under suitable additional conditions on {uo,vo), 

(5.5) limsup—^ < c*, 

where c* is the minimal speed of the traveling waves to (jl.5p in dimension one given in |8j . More 
precisely, by Theorem 2.1 of [8], when (j4.ip holds, there exists c* > 0, depending continuously on 
the parameters in (jl.5p . such that (jl.5p with = 1 has a solution of the form (jl.6p when c > c*, 
and there is no such solution when c < c=k. (We remark that by a standard change of variables, 
our general form (jl.5p can be reduced to the special form considered in [8].) 



Theorem 5.3. Suppose that (14. 1|) holds, h^o = 00, and additionally 

(5.6) uo(r) <^ in [0,ho], vo{r) > in [0,oo), liminf^^oo t'o(?') > 

Then (j5.5p holds. 
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Proof: In (jl.5p . we replace ai by ai = ai + ei and replace 02 by 02 = (12 — €2, where ei, £2 
are small positive constants. We then denote by c* the minimal speed of the traveling waves to 
the modified (jl.Sp . and by {U{x),V{x)) the corresponding solution in (jl.6p . We thus have 

' U" + cM' + U{ai-biU -ciV) = 0, U'<0, x G M^, 

< V" + c^V' + V{d2-b2U -C2V) = 0, V'>0, xGR\ 
^ U{-oo) = |, y(-oo) = 0, Ui+00) = 0, V{+oo) = |. 

Let ct(x) be a smooth function satisfying 

a{x) G [0, 1], f7(x) = for X < -1, a{x) = 1 for x > 0, a'{x) > 0. 

For xo > and L > to be determined later, we define 



U{x 

Clearly U{x) < U{x) and 



r ;7(x)-a(2^^)C/(xo), x<xo, 
1 0, X > Xo- 



U'{x) < 0, X G M\ 

lJ{x) = U{x) -U{xq), xG[xo-1,xo], 
ij{x) = U{x), x<xq-1-L. 

We now define 

m = c,t + ^o, u*{t,r) = U{r- C{t) + xq), v,{t,r) = V{r - ^{t) + xq). 

We will show that by choosing xo,L,S^q properly, we have h{t) < ^(t) for all t > 0. Clearly this 
implies 

lim sup < lim = . 

t^oo t t->-oo t 

Letting ei, £2 — )• 0, due to the continuous dependence of c* on the parameters, we deduce 

h(t) 

lim sup < c* , 

as desired. 

It remains to show h{t) < ^(t) for t > 0. This will be shown by checking that (n*,^^,,.^) can 
serve as an upper solution to (jl.ip . Clearly 

u*{t,C{t)) = U{xo) = > 0. 

Moreover, 

= /xf/'(xo) = fiU'ixo) < = c'{t) 

provided xq > is chosen large enough, since U'{x) — )• as x — )• +00. 
By direct calculations we obtain 

— diAu* 

= -cM'{r - m + xo) - di (u"{r - ^{t) + xq) + ^^^U'{r - ^{t) + xq)^ 
> -cM' - diU" 

c* — + di 22- j 

a" 



Uiai - biU - ciV) + L-^Uixo) i^c^ + di ^ 
U{ai - biU - ciV) + eiU + L-^U{xq) {c^a' + c^i^) • 
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For r G (0, ^(t)], we have r — ^{t)+XQ < Xq and therefore U{r — ^,{t) + XQ) > U (xq). On the other 



hand, by choosing L suitably large we can guarantee that L ^ 
L chosen this way, we have 



c*(t' + d 



1 — 



< ei. Therefore, with 



eiU + L-^U{xo) i^c,a' + di^j > for r G [0,C{t)]. 

Moreover, we notice that U{x) = U{x) for x < xq — 1 — L, and by enlarging xq if necessary, we can 
make {U{xo — 1 — L), V{xo — 1 — L)) very close to (0, ^). It follows that for x £ [xq — 1 — L, xq], 

(f7(x), V{x)) is very close to (0, ^). Hence, due to (|4.ip . for every fixed x G [xq — 1 — L,xq], the 
function 

r]{s) := s[ai — bis — ciV{x)] is increasing in [0,U{x)]. 
We may now use U{x) > U (x) to deduce that, for x £ [xq — 1 — L, xq], 

U{x)[ai - biU{x) - ciV{x)] > U(x)[ai - biU{x) - ciV{x)]. 
We thus obtain, by combining the above inequalities, 

ul — diAti* > U{ai — biU — ciV) = u*{ai — biu* — civ^) 
for r £ (0, Cit)] and t > 0. 



We also have 



- 1 

{v,)t - d2Av, = -c,V' - d2V" - d2 V 



< -c^V - d2V" 

= V{d2 - b2U - C2V) 

< Vid2 - b2U - C2V) 

< V^{a2 - 62 It* - C2f*) 

for r > and t > 0. 

For r € [0, /iq], u*{0, r) = U{r — + xq) > U{hQ — ^0 + xq). Since U{x) — )• as x — )■ —00, we 
may choose .^0 > ^0 large enough such that U {ho — ^0 + xq) > ^ and hence, with such a .^Oi 

u*{0,r) > u{0,r), Vr G [0,ho]. 

By our assumption on V(){r), there exists Xi > such that VQ{r) > ^ for r > Xi. Let 
(To := infj,g[o^Xi] vo{r). Then uo > and we can find X2 > such that V{x) < ctq for x < —X2. 
Therefore, for r £ [0,Xi], 

v^O, r) = V{r - ^0 + xq) < V{Xi - + xq) < ctq < vo{r) 
if we enlarge .^0 further to guarantee Xi — S^q + xq < —X2. Since w=„(0, r) = V{r — ^0 + xq) < ^ 
for every r, and VQ{r) > |j for r > Xi, we thus have 

i'*(0,r) < vo{r) Vr > 0. 

We are now in a position to apply Lemma [2.61 to deduce h{t) < ^{t), u{t,r) < u*{t,r) and 
v{t,r) > 'L'*(t,r), except that the conditions drU*{t,0) = 5r'L'*(t, 0) = are not satisfied. Instead, 
we have 

drU*{t,0) = U' {-C^t - + Xq) < 0, 5,.W*(t,0) = V {-C^t - + Xq) > 0. 

This case is covered by Remark \2.7\ and hence we can still apply the comparison principle to 
obtain the stated inequalities. □ 

Remark 5.4. It is easily seen from the above proof that if instead of (|5.6|) . we assume that there 
exists T > such that u{T,r) and v{T,r) satisfy the conditions in ()5.6p . then the conclusion of 
Theorem 15.31 remains valid. 
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